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A universal quantum work relation is proved for isolated time-dependent Hamiltonian systems 
in a magnetic field as the consequence of microreversibility. This relation involves a functional of 
an arbitrary observable. The quantum Jarzynski equality is recovered in the case this observable 
vanishes. The Green-Kubo formula and the Casimir-Onsager reciprocity relations are deduced 
thereof in the linear response regime. 

Nonequilibrium work relations have recently attracted much interest Q . They provide relations for the work 
dissipated in time-dependent driven systems, independently of the form of the driving. They are of great interest to 
evaluate free energies under general nonequilibrium conditions and they provide new methods to study nanosystems. 
In the nanoscopic world, the extension of these classical relations to quantum systems is of particular importance and 
different approaches have been proposed. 

A first scheme was introduced by Kurchan [s^ . In this framework, a measurement of the system state is performed 
at the initial time. In the sequel, the system is perturbed by a time-dependent Hamiltonian before performing another 
measurement at the final time. The random work performed on the system is associated with the energy difference 
between the final and initial eigenstates. This setup leads to the quantum extension of Jarzynski equality and Crooks 
fluctuation theorem 0, 0, Q • Another possibility is to introduce a quantum work operator which measures the 
energy difference Q, in which cases quantum corrections to the fluctuation theorem must be taken into account. On 
the other hand, quantum fluctuation theorems have been obtained in suitable limits where the dynamics admits a 
Markovian description, allowing in particular the applications to nonequilibrium steady states [ol. [lOl [ill, [l^. fisl. [T^ . 
Yet, the connection between the quantum work relations and response theory is still an open question even in the 
linear regime. 

The purpose of the present paper is to derive a new type of work relations which involves a functional of an 
arbitrary observable. This generating functional can be related to another functional but averaged over the time- 
reversed process. This new work relation turns out to be of great generality since we can recover known results such 
as Jarzynski equality as special cases. Furthermore, this universal work relation allows us to formulate the response 
theory, to derive the quantum linear resp onse functions, the quantum Green-Kubo relations [Tsl. [T6|. as well as the 
Casimir-Onsager reciprocity relations [l3, [3 in the regime close to the thermodynamic equilibrium. 

Functional symmetry relations. We suppose that the system is described by a Hamiltonian operator H{t]B) which 
depends on the time t and the magnetic field B. The time-reversal operator O is an antilinear operator such that 
= / and which has the effect of changing the sign of all odd parameters such as magnetic fields: 

eH{t;B)e^ H{t;~B). (1) 

We first introduce the forward process. The system is initially in thermal equilibrium at the inverse temperature 
f3 — 1/k^T. The initial state of the system is described by the canonical density matrix 

-/3H(0:e) 

where the partition function is given in terms of the corresponding free energy F(0) by Z(0) = tr e^'^^'^*'-^) = e^^^("^ 
Starting from this equilibrium situation at the initial time t = 0, the system evolves until some final time t ^ T under 
the Hamiltonian dynamics. The corresponding forward time evolution is defined as 

ih^UF{t;B) = H{t-B)UF{t-B), (3) 
with the initial condition U-p{Q]B) = I [l^. In the Heisenberg representation, the observables evolve according to 

AF{t)^Ul{t)AUF{t) (4) 
which also concerns the time-dependent Hamiltonian 

HF{t)^Ul{t)H{t;B)UF{t). (5) 
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The average of an observable is thus obtained from 

(AfW) -trp(0)AF(i). (6) 

We note that the dependence on the magnetic field is implicit in these expressions. 

The backward process is introduced similarly but in the magnetic field reversed. The system is perturbed according 
to the time-reversed protocol H(T — t; —B), starting at the initial time t — from the density matrix 

-I3H{T:-B) 

PiT)-^^, (7) 

where the free energy F{T) is given in terms of the partition hmction according to Z{T) — tre"'^^^"^'"^-' = e~^^'^'^\ 
The system ends at time t — T with the Hamiltonian -ff(0; —B). The evolution operator of the backward process is 
defined as 

ih-UR{t;B)=HiT-t;B)Unit;B), (8) 

with the initial condition L/r(0;B) = / [l^, and is related to the one of the forward process by the following 

Lemma: The forward and backward time evolution operators are related to each other according to 

QUf{T ~ t; B)UI{T- B)Q = Un{t; -B) , (9) 
where t is an arbitrary time < t < T . 

This lemma is proved by first substituting T — t for t in Eq. ^ to get 

- ih^UpiT -t;B)^ H{T - t; B)Uf{T -t;B). (10) 
Multiplying this equation by Ul{T; B)Q from the right and by G from the left, we find 

ih^QUFiT ~t-B)Ul{T-B)Q = H{T - t; -B)eUFiT ~ t; B)Ul{T; B)Q , (11) 

where we used the antilinearity Oi — —iQ of the time- reversal operator and its further property ([T]). This shows that 
the expression QUf{T — t; B)Ul{T; B)Q obeys the same evolution equation ^ as U-R,{t; —B). Since they also satisfy 
the same initial condition, eUpiT; B)Ul{T; B)Q = Ur{Q; -B) = /, we have proven Eq. dH). QED. 

With this lemma, we can now demonstrate the 

Theorem: Let us consider an arbitrary time-independent observable A with a definite parity under time reversal: 
QAQ = caA, with = ±1- It satisfies the following functional relation: 

j;^ dtX{t)AF(t) -I3Hf{T) I3H{0)\ ^ P-/3AF / f;f dtX{T^t)An{t) 



. , (12) 

/F,B \ /R,-6 

where X{t) is an arbitrary function, while the subscripts F and R stand for the forward or backward protocol, 
respectively. AF = F(T) — F(0) is the difference of the free energies of the initial equilibrium states (Qj and 0) of 
the backward and forward processes. 

In order to prove Eq. ()12p . we first consider the quantity Aplt), which can be written as 

AF{t) - ul{t)AUF{t) - uliT)UFiT)ul{t) AUF{t)ul{T)UF{T) = eAUl{T) Ar{T - t) eUF{T) , (13) 

where we have inserted the identity uI{T)Uf{T) = / to go at the second equality. At the third equality, we inserted 
= I between the evolution operators and we used QAQ = caA along with Eq. The connection is thus estab- 
lished with the backward process. Integrating over time with an arbitrary function X{t) and taking the exponential 
of both sides, the previous expression becomes 

exp dtX{t)AFit)^=Ul{T)Qexp(^eAj^ dt X(T ~ t) A^it)^ Q Uf{T) , (14) 
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after the change of integration variables t — s- T — t in the right-hand side. 
Starting from the left-hand side of Eq. we get 

trp(0)exp dtX{t)AF{t)^ exp[-/?i/F(T)] exp[/3i/(0)] 

= ^ tr exp (^SA ^ dt A(T - t)An{t)^ 6 exp[-/3£f (T; B)]e 



Z{T) 
Z{0) 



trexpfe^^ A(T ~ i)^R(t) j p(T) = e"^^^ /exp ^ A(r - t)AR,(i) j \ ■ (15) 



R,-8 



We used the invariance of the trace over cyclic permutations as well as the exponential of Eq. (fT^)) at the first equality. 
In the second equality, we introduced the equilibrium density matrix ([7]) which is precisely the initial condition of 
the backward process. To obtain the last equality, we used that the partition functions have been expressed in terms 
of the corresponding free energies. This completes the proof of the theorem. QED. 

We notice that related results have previously been considered in the restricted case where there is no change in 
free energy AF = [20l . [2l| . The present theorem allows us to recover in particular the quantum Jarzynski equality 
as a special case of Eq. ([T^ if A = 0: 



/f.b 



The factor inside the bracket can indeed be interpreted in the quantum setting in terms of the work performed 
on the system during the forward process [1, 0, [a, Q in spite of the non-commutativity of the energy operators 
H-p{T) and i?(0) and thanks to the protocol with von Neumann quantum measurements of the energy at the initial 
and final times. It is only in the classical limit that both energies commute and the classical work can be formed as 
Wci = [H-p{T) — i7(0)]j,j. In this case, both exponentials in the left-hand side of the relation p2|) becomes exp(— /3Wci) 
which is the classical version of this relation. 

Response theory. We can obtain different correlation functions by taking functional derivatives of the relation (jl2p 
with respect to the arbitrary function A(<). In this way, we can obtain the expression of linear response theory from 
the generalized symmetry relation (jl2p . For this purpose, we consider a perturbation of the form 

H{t) = Ha-X{t)B, (17) 

where the perturbation X(t) is such that X(i) = for t < and X{t) = for T < t. The observable B is here 
arbitrary and should not be confused with the magnetic field B. In order to obtain the linear response of an observable 
A with respect to the perturbation —X{t)B, we take the functional derivative of Eq. with respect to A(r), around 
A = 0. This yields 

(AAT)e-^"^(^^e^"A = (Ar(0))j,__^ = {A)^^ _^ , (18) 

where we used that AF = since X{0) = X{T) — 0. Since the reversed process also starts at equilibrium, the 
average in the right-hand side is an equilibrium average, albeit with a reversed magnetic field. Nevertheless, we have 
that eA (^)oq ~ (^)cq B using time reversal. We now have to calculate the exponentials of the initial and final 
Hamiltonians. Since, in the Heisenberg representation, the total time derivative of the Hamiltonian equals its partial 
derivative, dHp/dt ~ {dH/dt)-p, we can write 

exp[-/3i/F(T)] = exp[-/3(i/o + E)] (19) 

with 

I dt(^\ =-( dtX{t)BF{t)^ [ dtX{t)BF{t), (20) 



(21) 



"'0 \ / F "'O "'0 

where the last equality follows from an integration by parts. We now use the expression 

r/3 



exp[/3(P -I- Q)] exp(-/3P) = 1 + / du exp[u(P + Q)] Q cxp(-uF) 

Jo 
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which can be proved by differentiating with respect to /3. To first order in Q, we may neglect Q in the last exponential 
function, exp[u(P + Q)]. Taking P — —Hq and Q — —E and developing to first order in X, we get 

g-/3ifp(T)g/3ifo = 1 _ / dtX{t) due-"^«S(t)c"^° +0(X2) = f - / dtX{t) duB{t + ihu) + OiX'^), 
Jo Jo Jo Jo 

where B{t) = eKp(iHot / h) B exp(—iHot/ h) since, at first order in the driving force, the time evolution proceeds under 
the unperturbed Hamiltonian Hq. Inserting this expansion into Eq. (|f 8p and after some manipulations using the time 
invariance of correlation function as well as the KMS-like property pA = A(ihj3)p [l^], we finally find 

{A^{T))^ = {A) + j dtX{T-t)^AB{t) + 0{X^) , (22) 
Jo 

with the response function 

(t>AB{t) = / du {B{-ihu)A{t)),q,B . (23) 
Jo 

Equations and P5)) are the well-known expressions of linear response theory in the canonical ensemble, also 
known as the Green-Kubo formula [H, • The Casimir-Onsager reciprocity relations for the conductivities EBl 
are obtained by taking A = J^/V and B = in terms of the current = J2n^ninfj. and the volume V, in 
which case the time-reversal symmetry implies (jj^i^lt; B) — (pi^fiit; —B) and a^^luj] B) — av^{uj; —B) for the tensor of 
conductivities cr^t^(w; B) = dte^^*(j)fj,,^{t; B). Higher-order terms in the expansion can be obtained as well. 

Conclusions. In this paper, we have obtained a universal quantum work relation which involves arbitrary 
observables at arbitrary times. This result relates an average over the forward process ponderated by the quantum 
analogue of the work to an average over the reversed process. By taking functional derivatives, we can obtain 
relations for arbitrary correlation functions, which are the consequence of microreversibility. In the simplest case, it 
can be used to recover the well-known Jarzynski equality. On the other hand, we can also straightforwardly derive 
from the universal relation the linear response theory of an arbitrary observable. In this regard, this relation unifies 
in a common framework the work relations and the response theory, thereby opening the possibility to obtain further 
general relations which are valid not only close to equilibrium but also in the far-from-equilibrium regime. 
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